Mass-imbalanced Fermi gases with spin-orbit coupling 
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We use the mean-field theory to analyze the ground-state phase diagrams of spin-orbit coupled mass- 
imbalanced Fermi gases throughout the BCS-BEC evolution, including both the population-balanced and - 
imbalanced systems. Our calculations show that the competition between the mass and population imbalance 
and the Rashba-type spin-orbit coupling (SOC) gives rise to very rich phase diagrams, involving normal, su- 
perfluid and phase separated regions. In addition, we find quantum phase transitions between the topologically 
trivial gapped superfluid and the nontrivial gapless superfluid phases, opening the way for the experimental 
observation of exotic phenomena with cold atom systems. 
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PACS numbers: 05.30.Fk, 03.75.Ss, 03.75.Hh 

Introduction. Following the great successes with single- 
species (mass-balanced) two-component Fermi gases Oil, 01, 
there has been increasing experimental interest in realiz- 
ing and studying two-species (mass-imbalanced) Fermi-Fermi 
mixtures over the last few years ^ . So far the most promi- 
nent candidate seems to be the ^Li-'^'^K mixtures, for which 
the experimental methods are currently being developed in 
several groups. For instance, ^Li-*°K mixtures have recently 
been trapped and interspecies Feshbach resonances have been 
identified, opening a new frontier in ultracold atom research 
to study exotic many -body phenomena. We also note that sev- 
eral other fermionic atoms including ^^^Yb |5] and ^^Sr Ig] 
are also cuiTently being investigated, allowing for future mix- 
ture experiments with various other species. Motivated partly 
by these experiments, and also due to the natural way of cre- 
ating superfluidity with mismatched Fermi surfaces, there has 
also been increasing theoretical interest in understanding and 
studying two-species mixtures ||7|,[al- 

In addition to these developments, following the recent 
realization of synthetic gauge fields with neutral bosonic 
atoms [9] and the more recent creation of spin-orbit coupled 
BECs 1 10], it is now possible to create and study spin-orbit 
coupled Fermi gases, by making use of similar experimen- 
tal methods |11]. Since this new technique allows for the 
realization of topologically nontrivial states in atomic sys- 
tems with possibly a broad interest in the physics commu- 
nity I.12.-.141. there has been increasing theoretical interest 
in studying the effects of SOC on the single-species two- 
component Fermi gases. For instance, it has been found for 
the population-balanced Fermi gases that the SOC increases 
the single-particle density of states, which in return favors 
the Cooper pairing so significantly that increasing the SOC, 
while the scattering length is held fixed, eventually induces 
a BCS-BEC crossover even for a weakly-interacting system 
when as -^ 0~ fllsMlTll . The increased density of states also 
has important effects on the population-imbalanced Fermi 
gases Ill8[[l9ll . for which we have recently found that the SOC 
and population imbalance are counteracting, and that this 
competition tends to stabilize the uniform superfluid phase 
against the phase separation. In addition, topological quan- 
tum phase transitions associated with the appearance of mo- 



mentum space regions with zero quasiparticle energies have 
been found, the signatures of which could be observed in the 
momentum distribution Ill8l Il9ll . 
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FIG. 1: (Color online) The ground-state phase diagrams are shown 
at unitarity, i.e. the scattering parameter is set to l/{kFas) = 0, as 
a function of population imbalance P — {Ni; ~ Ni)/N and mass 
ratio m^/mj_. Here the SOC parameter am+/kF is set to 0.05 in 
(a), 0.15 in (b), 0.3 in (c), and 0.45 in (d), where m+ is twice the 
reduced mass of ^ and J, fermions. We show normal (N), phase sepa- 
ration (PS), topologically trivial gapped superfluid (SF) and the topo- 
logically distinct gapless superfluid (GSF) phases. 

In this paper, we extend our recent work lUal to the case 
of mass- and population-imbalanced Fermi gases, and study 
the effects of spin-orbit coupling on the ground-state phase 
diagrams across a Feshbach resonance, i.e. throughout the 
BCS-BEC evolution. Since it is not possible to have a 
spin-orbit field that converts different species of atoms into 
each other, it may not be possible to create spin-orbit cou- 
pled mass-imbalanced systems with two-species mixtures. 
However, mass-imbalanced systems can be engineered with 
single-species two-component Fermi gases loaded into spin- 
dependent optical lattices in such a way that the components 



have different effective masses. 

Some of our main results are shown in Fig. [T] and they 
are as follows. The competition between the mass and pop- 
ulation imbalance and the SOC gives rise to very rich phase 
diagrams, involving normal, superfluid and phase separated 
regions, and quantum phase transitions between the topolog- 
ically trivial gapped superfluid and the nontrivial gapless su- 
perfluid phases. In addition, in sharp contrast to the no-SOC 
case where only the gapless superfluid phase can support pop- 
ulation imbalance, both the gapless and gapped superfluid 
phases can support population imbalance in the presence of 
a SOC for all mass ratios including the mass-balanced Fermi 
gases. Similarly, in again sharp contrast to the no-SOC case 
where only the gapped superfluid phase can support popula- 
tion balance, both the gapped and gapless superfluid phases 
can support population balance in the presence of a SOC when 
the mass difference becomes large enough. 

Mean-field theory. We obtain these results within the self- 
consistent mean-field approximation. In the absence of a SOC 
and at low temperatures, it is well-established that the mean- 
field theory is sufficient to describe the physics of Fermi gases 
both in the BCS and the BEC Umits, and that this theory also 
captures qualitatively the correct physics in the entire BCS- 
BEC evolution 1 1, 20]. Hoping that the mean-field formalism 
remains sufficient in the presence of a SOC, which is expected 
since the system is weakly interacting both in the BCS and 
molecular BEC limits, here we analyze the ground-state phase 
diagrams of the system as a function of SOC, scattering and 
population imbalance parameters as well as the mass ratio of 
the fermions. 

For this purpose, we use the mean-field Hamiltonian (in 
units of h = 1 — ks) 
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(flk.cr) creates (annihi- 



lates) a spin-cr fermion with momentum k. Here, C = 

(1/2) Ek.CT '?k,ff + I A|2/g is a constant, and fk,<T = ek.^ - f^a 
with ek,o- = k"^ /{2ma) the kinetic energy, /Xo- the chemi- 
cal potential and k = Jk'^ + k^ + k'^. In addition, 5k — 

a{ky — ikx) with strength a > is the Rashba-type SOC ll2lll . 
and A — .g(ak,ta_k,4.) is the mean-field order parameter, 
where (• • • ) is a thermal average and 5 > is the strength 
of the attractive particle-particle interaction which is assumed 
to be local. 

The corresponding mean-field thermodynamic potential 
can be written as 
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where T is the temperature, s — ±, and E^ ^ = Ck + + ^ 



k. 



trum MM. Here, A^ = ^JClA^i+ + |AP) + l^kP^^ + 

and fk,s = ek,s - Ms, where ek,s = (ek,t + ■sek,i)/2 = 
k'^/{2ms), rUs — 27Ti-|-?Ti_[,/(77ij_ + snif) and /i<, = (/i^ + 
s/i4.)/2. Following the usual procedure, i.e. dft/d\A\ = for 
the order parameter and A^-|- + sN^ = —dil/dfis for the num- 
ber equations, we obtain the self-consistency equations 1 18] 
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Here, we eliminated the theoretical parameter g in favor of the 
experimentally relevant s-wave scattering length a^ via the re- 
lation, 1/5 = —raj^V/{ATTas) + X^k l/(2ek.+), where m+ is 
twice the reduced mass of f and J, fermions and V is the vol- 
ume. The derivatives of the quasiparticle energies are given 
bya£;k,s/9|A|2 == (l + s^g _/ylk)/(2£;k,s) for the order pa- 
rameter, aSk,s/9M+ = -[l + s(Ck,- + l^kP)Mk]ek,-/-Bk,s 
for the average chemical potential and 9i?k.s/9/i- = — [1 + 
*(^k + + I AP)Mk]Ck,-/-Ek,s for the half of the chemical po- 
tential difference. 

We checked the stability of the mean-field solutions for 
the uniform superfluid phase using the curvature criterion 13], 
which says that the curvature of Vt with respect to | A|, i.e. 




I^kp + 2sy4k gives the quasiparticle excitation spec- 



(5) 



needs to be positive. When the curvature d'^il/d\/S.\'^ is neg- 
ative, the uniform mean-field solution does not correspond to 
a minimum of 57, and a nonuniform superfluid phase, e.g. 
a phase separation, is favored. It is known that the curva- 
ture criterion correctly discards the unstable solutions, but the 
metastable solutions may still survive. This may cause only 
minor quantitative changes in the first order phase transition 
boundaries as shown for the mass-balanced Fermi gases 111911 . 
Having set up the mean-field formalism, we are now ready to 
discuss the competition between normal fluidity, uniform su- 
perfluidity and phase separation across a Feshbach resonance. 
Ground-state phase diagrams. There are typically three 
phases in our phase diagrams. While the normal (N) phase is 
characterized by A = 0, the uniform superfluid and nonuni- 
form superfluid, e.g. phase separation (PS), are character- 
ized by (92f7/a|A|2 > and d'^n/d\/^\^ < 0, respectively, 
when A 7^ 0. Furthermore, in addition to the topologically 
trivial gapped superfluid (SF) phase, the gapless superfluid 
(GSF) phase can also be distinguished by the momentum- 
space topology of its excitations. Depending on the number 
of zeros of _Ek.s (zero energy points in k-space), there are 
two topologically distinct GSF phases: GSF(I) where E^s has 
two, and GSF(II) where E^g has four zeros. For the Rashba- 
type SOC, the zeros occur when k± — Q and at real kz mo- 



menta, k'^ g = -B+ + sJB'^ — 4?n-|-m^|Ap, provided that 

|A|2 < |B_|V(4mtm;) for B+ > 0, and jAp < -fi-^fii 
for B+ < 0. Here, Bg = m^fi^ + smifii- The topologically 
trivial SF phase corresponds to the case where both -Ek,+ and 
i?k,- have no zeros and are always gapped. From this analy- 
sis, it follows that the conditions |Ap = |i?_|^/(4m-|-m4,) and 
Ap = —jj^iJi determine the phase boundaries between the 
SF, GSF(I) and GSF(II) regions, and that these three phases 
meet at a tricritical point determined by B+ = 1 18]. 
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FIG. 2: (Color online) The ground-state phase diagrams are shown 
as a function of m^/mi and a at l/(fcFas) ~ 0, where P is set to 
0.25 in (a), 0.5 in (b), -0.25 in (c), and -0.5 in (d). The labels ai-e 
described both in Fig.[T]and in the text. 

In our numerical calculations, we use an effective Fermi 
momentum kp and the corresponding Fermi energy ep — 
kp/{2m+) as our length and energy scales, where kp is 
determined by fixing the total number of fermions N = 
N-f + Ni = A:|y/(37r^). In addition, we choose t ii) to 
label lighter (heavier) fermions such that lighter (heavier) 
fermions are in excess when the population imbalance param- 
eter P — {N-f — Ni)/N is positive (negative). Therefore, 
since we choose — 1 < P < 1, it is possible to span all pos- 
sible mass-imbalanced Fermi gases by considering the mass 
ratios < m-|-/m^ < 1. 

(I) Generel phase diagrams. In Figs. [T] and |2] the ground- 
state phase diagrams are shown as a function of P and m-^/mi 
for fixed a values, and as a function of a and m^/ra^ for 
fixed P values, respectively, at unitarity, i.e. \/{kpas) = 0. 
The dashed blue and dotted green lines correspond to | A|^ — 
— /i^/ij. and |A|^ = |i?_p/(4TO-t-77T,|), respectively. Since our 
classification of distinct topological phases applies only to the 
uniform superfluid region, the dashed and dotted lines shown 
within the PS regions are solely for illustration purposes. 

We find that the phase diagrams are symmetric around 
P — Q for mass-balanced systems, and that this symmetry is 
gradually broken as the mass difference increases. Due to this 
asymmetry, while the N, SF and GSF(II) phases occupy much 
larger regions when heavier fermions are in excess, the PS and 



GSF(I) phases occupy larger regions when lighter fermions 
are in excess. It is clearly shown in these figures that increas- 
ing a gradually stabilizes the SF and GSF phases against the 
N and PS for all mass ratios, and mostly the SF and GSF(I) 
phases remain in the phase diagrams for very large a. This is 
mainly a consequence of increased single-particle density of 
states due to the SOC as mentioned in the introduction. 

In sharp contrast to the a — Q case where only the gap- 
less GSF phase can support population imbalance, one of the 
intriguing effects of the SOC is that both the gapless GSF 
and gapped SF phases can support population imbalance when 
a 7^ 0. As shown in Figs. [T] and |2] this happens for all mass 
ratios including the mass-balanced systems |18, 19], but the 
gapped SF phase with population imbalance occupies a much 
larger region when the heavier fermions are in excess. Sim- 
ilarly, in again sharp contrast to the a = case where only 
the gapped SF phase can support population balance, another 
intriguing effect of the SOC is that both the gapped SF and 
gapless GSF phases can support population balance for large 
enough mass differences when a ^ 0. At unitarity, this occurs 
for all a values when m^ < 0.5to|, with the largest effect at 
around a ~ 0.55kp/m^. 
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FIG. 3: (Color online) The ground-state phase diagrams of a Fermi 
gas with m-f — 0.15771^ are shown as a function of P and l/{kFas), 
where amj^/kp is set to 0.05 in (a), 0.15 in (b), 0.3 in (c), and 0.45 
in (d). The labels are described both in Fig.[T]and in the text. 

(//) Fermi gas with m^ = 0.15?7i^. Having discussed the 
general structure of the phase diagrams as a function of mass 
ratio, here we fix nif = 0.15m^ and present a thorough dis- 
cussion of the resultant phase diagrams in the entire BCS-BEC 
evolution. This particular mass ratio corresponds to that of 
two-species ^Li-^'^K mixtures which are currently being de- 
veloped in several groups, making them the most prominent 
candidate for realizing a mass-imbalanced Fermi gas yfl- 

Our results are shown in Figs.[3]and|4l where the ground- 
state phase diagrams are shown as a function of P and 
l/{kpas) for fixed a values, and as a function of P and a for 
fixed l/{kpas) values, respectively. In addition to the general 
findings discussed above, we find that the gapped SF phase 



with population imbalance gradually disappears towards the 
BEC limit giving its way to the gapless GSF(I) phase. This 
is quite intuitive since, in this limit, the Fermi gas is expected 
to be well-described by a much simpler Bose-Fermi descrip- 
tion of paired fermions (molecular bosons) and unpaired (ex- 
cess) fermions, consistent with the k-space topology of the 
GSF(I) phase [7]. In addition, we find that both the gapped 
SF and gapless GSF(II) phases occupy much larger regions 
for all l/{kpas) values when heavier fermions are in excess. 
When m^ = O.lBm^, we finally note that the gapless GSF(I) 
phase with population balance can be observed at unitarity 
for moderate a values, through for instance measurements of 
the momentum distribution and single-particle spectral func- 
tion lal. 




0.25 0.5 0.75 1 0.25 0.5 0.75 1 
am^jkp am^/kf 

FIG. 4: (Color online) The ground-state phase diagrams of a Fermi 
gas with m-|- = O.lSmj, are shown as a function of P and a, where 
l/ikpas) is set to —1 in (a), —0.5 in (b), in (c), and 0.5 in (d). The 
labels are described both in Fig.[T]and in the text. 



Conclusions. To summarize, we analyzed the effects of 
SOC on the ground-state phase diagrams of mass-imbalanced 
Fermi gases throughout the BCS-BEC evolution. One way to 
engineer such a system is to load a two-component Fermi gas 
into a spin-dependent optical lattice such that the components 
have different effective masses. We showed that the com- 
petition between the mass and population imbalance and the 
SOC gives rise to very rich phase diagrams, involving normal, 
superfluid and phase separated regions, and quantum phase 
transitions between the topologically trivial gapped superfluid 
and the nontrivial gapless superfluid phases. According to our 
calculations, one of the intriguing effects of the SOC is that, 
in sharp contrast to the no-SOC case where only the gapless 
superfluid phase can support population imbalance, both the 
gapless and gapped superfluid phases can support population 
imbalance in the presence of a SOC. This occurs for all mass 
ratios including the mass-balanced systems. Another intrigu- 
ing effect of the SOC is that, in again sharp contrast to the no- 
SOC case where only the gapped superfluid phase can support 
population balance, both the gapped and gapless superfluid 
phases can support population balance in the presence of a 
SOC when the mass difference becomes large enough. 
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